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Abstract 

We examine the space-time evolution of (pre-)hadron production within the Lund 
string fragmentation model. The complete four-dimensional information of the string 
breaking vertices and the meeting points of the prehadron constituents are extracted 
for each single event in Monte Carlo simulations using the jETSET-part of Pythia. 
We discuss the implication on the deep inelastic lepton scattering experiments at 
HERMES as well as on observables in ultra-relativistic heavy ion collisions at RHIC, 
using Pythia also for modeling the hard part of the interaction. 
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The fragmentation of high energy quarks into color-neutral hadrons is a sub- 
ject of high current interest. So far one cannot describe the complete space- 
time picture of this process within QCD as the underlying theory. The reason 
is that the formation of hadrons involves momentum scales of only a few 
hundred MeV which rules out a purely perturbative approach. Furthermore, 
the complexity of the fragmentation process goes far beyond what can be 
addressed by current non-perturbative lattice-QCD calculations. 

Estimates based on perturbative QCD (jl|) suggest, that the formation time 
of a hadron, i.e. the time needed to build up the hadronic wave function, is 
of the order of Ifm/c in the hadron's restframe. At large enough energies 
these eigentimes correspond to formation lengths in the lab frame which are 
comparable with nuclear dimensions. If the fragmentation process is modi- 
fied by the surrounding nuclear medium, an experimental comparison of high 
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energy hadron production on various nuclear targets with different radii can 
give insight into the space-time picture of hadronization. With that respect 
the nucleus can be considered as a micro-laboratory for fragmentation studies. 



In the recent past the HERMES collaboration at DESY has carried out an ex- 
tensive experimental study of hadron attenuation (0) in deep inelastic lepton- 
nucleus scattering (DIS). Depending on the evolution of the fragmentation 
process in space-time the observed signals might either be dominated by (pre- 
)hadronic final state interactions (FSI) (0; 0; S 0) or by a partonic energy 
loss prior to fragmentation (0). A detailed investigation of this subject is also 
essential for the interpretation of ongoing experiments at hadron colliders 
like RHIC or the future LHC which among other topics intend to find ex- 
perimental evidence for a new QCD-phase, the Quark-Gluon-Plasma (QGP). 
(Pre-)hadronic FSI could contribute to the suppression of high-p^ hadrons 
and thereby modify jet-like signatures that are_still considered clean proofs 
for the creation of a deconfined QGP phase 



The Monte Carlo event generator Pythia ffiol ) has become one of the stan- 
dard tools toperform calculations in the kinematic regimes that we have just 



mentioned (la la; I 111 ). The fragmentation in Pythia is based on the Lund 
string model (|l2t). In former versions of Pythia the routines handling the 
string fragmentation were bundled in a separate package called Jetset. We 
keep this (old fashioned) nomenclature throughout this work in order to dis- 
tinguish between the part of the event generator that creates the hadronic 
strings (e.g. Pythia, Fritiof (flih etc.) and the part that is responsible for 
the string fragmentation (Jetset). A detailed description of the underlying 
model of Jetset can be found in Ref. (jl4j). In Jetset the eigentimes of 
the string breaking vertices are determined during the fragmentation. How- 
ever, the real essential information, namely the four- dimensional space-time 
coordinates of the string breaking, is not reported. 



In this work we present a detailed quantitative analysis of the four-dimensional 
space-time picture of hadron production with Pythia/ Jetset in soft and 
hard interactions at HERMES and RHIC energies. We show for the first time 
the exact spatiotemporal distribution of (pre-) hadron production and forma- 
tion points in realistic ep and pp collisions at HERMES and RHIC. Thereby 
we take into account complications like resolved photon interactions, gluon 
radiation, primordial transverse momenta of quarks, transverse momentum 
generation in the string fragmentation, cluster decay, finite mass effects, quark 
flavor dependence, diquark production, etc. Although our results should not be 
overstressed since one applies a semi-classical picture to a quantum mechan- 
ical problem, our approach makes it possible to assign the four-dimensional 
production and formation point to each single hadron in each single scattering 
event. This sets the foundation for consistently using the Pythia model for 
both the determination of hadron momentum spectra - where the model has 
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been successfully tested - and the description of the space-time evolution of 
hadron production. Both information can then be used in a consistent way 
as input for future transport simulations such as the BUU (5) or the HSD 
model (P) to describe nuclear reactions at HERMES and RHIC energies and 
to determine the contribution of the (pre-)hadronic FSI to the experimen- 
tally observed hadron attenuation. Since the application of our findings is not 
restricted to one single nuclear reaction model, we exclude any specific mod- 
eling of the (pre-)hadronic FSI at this stage and merely concentrate on the 
space-time evolution of the (vacuum) fragmentation process in this work. 

Although we use Pythia to produce the initial string configurations we em- 
phasize, that this is no restriction concerning the extraction of the four- 
dimensional fragmentation coordinates, since the fragmentation routines are 
independent of the string-building routines. It is therefore "ab initio" straight- 
forward, to transfer our method to other event generators, supposed they use 
Jetset for fragmentation. 



In principle, one could distinguish two classes of 'strings': strings that are made 
up of a quark and an antiquark only and those with additional (hard) gluons 
located between the string ends as transversal excitations. (We will not touch 
strings build up entirely of gluons in this work. Here the mechanisms at work 
are far beyond being transparent. On the other hand, they are very rarely 
produced at the HERMES energies and also play a minor role interpreting 
current RHIC findings.) However, this distinction is purely artificial, because 
the latter reduces to the first case in the limit of vanishing transverse gluon 
momenta We point out that Pythia also allows for diquarks at the 

string ends as well as the creation of diquark-antidiquark pairs in the string 
fragmentation. This possibility is accounted for in our computer code but 
omitted for the sake of better readability in the following text. 



Jetset includes two different hadronization methods for an initial partonic 
system: 'string fragmentation' and 'cluster decay' where only the first one is 
based on the actual Lund fragmentation model. The latter method is con- 
structed in such a way that it reproduces the rapidity distributions of the 
Lund fragmentation. Which of these two methods is used for hadronization 
depends mainly on the invariant mass of the decaying system. 

We begin with a review of the fragmentation of an one-dimensional string, 
i.e. a string that consists of a quark and an antiquark string end and no gluonic 
excitation. Note that, most authors that discuss 'string fragmentation' restrict 
themselves to this one type of all possibihties. It can easily be illustrated and 
analytic expressions can be found, e.g. in Refs. (0; [l^; lish . We thus explain 
the necessary formalism by means of this example (cf. Fig. 1). We stress, that 
all momenta and coordinates in Fig. 1 are four-vectors and not light-cone 
variables. 
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Fig. 1. A sketch of the fragmentation of an one- dimensional qq-String (in the rest 
frame of the string). The Px denote momenta, while Vi indicate string breaking 
vertices (production points) and the Fi the yoyo- formation points of the hadrons. 
While in principle all vectors are four- dimensional, in this sketch the horizontal 
axis is the x-axis and the vertical coordinate is the time t. Please note: not lightcone 
formalism, but four-vectors. 

In the following, we set the string tension to k = 1 GeV/ fm. This allows us 
to omit K, in all formulae and to identify momentum vectors with space-time 
coordinates, i.e. 

= x+ , Pg = X- (1) 

where Pq (Pq) and (x^) denote the initial four- momentum and turning- 
point of the (anti-)quark respectively. The four-momenta Pi of the final hadrons 
are given by 

P^=Cx^+^iX~ (2) 

with the constraints < < 1 and Y^^t = 1- For each hadron and 
^~ are obviously related via the (tranverse) mass rrii of the produced hadron. 
In Eq. (2) we have neglected any transverse momenta created in the string 
breakings for simplicity. The following formulae are not affected by this simpli- 
fication. The space-time positions of the string breakings (production points) 
are 

K = xtx^ + x7x- (3) 

with = 1 — I]}=i^/) X7 = + J2)=i^j^ and the points where the world- 
lines of the constituents of hadron i meet for the first time - usually called 
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the 'yoyo-formation point' of hadron i - are given as 



xtx~^ + Xi+iX . (4) 



One should keep in mind, that the latter only makes sense if one neglects any 
transverse momentum of the gg-pairs created at the string breaking vertices. In 
general, a quark from one vertex - having some transverse momentum - does 
not cross the trajectory of an antiquark that stems from another vertex and 
has some other transverse momentum. We point out that the yoyo-formation 
points do not necessarily correspond to the real formation points of the hadron, 
i.e. the space-time point where the full hadronic wavefunction has built up. 

Obviously, both the production and formation points in space-time can be 
reconstructed from the initial quark and antiquark momenta Pg and Pq if one 
extracts the complete set of from Jetset. (Please recall: All the transver- 
sal momentum components are only necessary to reconstruct the hadron mo- 
menta, not the fragmentation points.) 

We now concentrate on the implementation of the string fragmentation in 
Jetset. Like all other Lund-implementations it suffers from the fact that in 
the derivation of the fragmentation function one assumes that the decaying 



string has infinite mass If one starts on one side of the string (e.g. the 

quark string end) and continues the fragmentation to the other side (i.e. the 
antiquark string end), one will run into the problem of finding a solution for 
the last hadron that both conserves energy and momentum and fulfills the 
mass-shell condition for the hadron. In order to avoid this problem, Jetset 
randomly performs fragmentations at one or the other side of the decaying 
string, until the invariant mass of the remaining string system (somewhere 
in the middle of the string) falls below some threshold value. Then Jetset 
switches to a mode in which the parameters of the two final hadrons are chosen 
simultaneously in such a way that energy and momentum are conserved and 
that the rapidity distribution of the hadrons looks like in the original Lund 
fragmentation process (jiot llGh . This method is very similar to the one used 
for the 'cluster decay' discussed below. 

So far, the production and formation points could be determined by Eqs. (3) 
and (4) respectively. However, strings may have transverse excitations, which 
are represented by additional gluons between the string ends. These gluons 
impose some transverse momentum on the string and lead to a complex move- 
ment of the string in space-time which is known as the 'dance of the butterfly' 
(0). 

The additional gluons split the string into separate regions. For a configuration 
that consists of a quark, n — 2 gluons and an antiquark with related four- 
momenta q{pi), g{p2), ■■■ , g{Pn-i), QiPn) the initial string contains n — 1 
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pieces. In four-momentum space the i-th string piece is spanned by the two 
four- momenta p^-j and where we have defined p^-, = pi, p^-, = P2/2, 

Pp) = P2/2, p^2) = P3/2, . . . , = Pn-i/2, = pn in analogy to the 

simple gg-string (Eq. (1)). The factors 1/2 arise because the gluons share their 
energy between the two adjacent string pieces. 



Due to the complex movement of the string new regions appear as time goes by. 
Each region can be thought of being spanned by one p^-j and another four- 
vector with j and k not necessarily adjacent. It is straightforward, but tedious 
and not illuminating, to extend Eqs. (2) and (4) to the case of multiple string 
pieces ()l4^. Therefore we only comment on the major complications: Since 
now the string(-part) ends have initial "transversal" momentum components, 
transverse momenta of the quarks originating from the string breakings must 
be considered explicitly. Hadrons built up from partons that are created in 
different string regions cause severe problems in the string fragmentation. 
Simplifications that are made in the Jetset code to cure this problems result 
in faulty four-dimensional coordinates of the fragmentation points (|lj). 



As pointed out not every string configuration in Jetset decays according to 
the Lund fragmentation scheme. Depending on the invariant mass of the sys- 
tem it may also fragment according to a cluster decay algorithm where the 
final state may only consist of one or two hadrons. In case of the decay into one 
hadron, the strings's four-momentum is mapped onto the four-momentum of 
the final hadron. Since the hadron mass is more or less fixed, left over compo- 
nents are transfered to a nearby string in form of a gluon with the correspond- 
ing four-momentum. If the cluster decays into two hadrons, the same algorithm 
as for the two last hadrons in the usual string fragmentation is applied. For de- 
tails we refer the reader to Ref. ^di) and references therein. Production vertices 
of hadrons stemming from cluster decays are assumed to be at the interaction 
point. In case of a cluster decay into two hadrons, the four-dimensional infor- 
mation of the second, i.e. the "middle" quark-antiquark production point, is 
reconstructed from energy- momentum constraints in analogy to Eq. (3). Also 
the formation points are determined from the four-momenta of the hadrons 
in analogy to Eq. (4). 

Since Jetset does not report any space-time vectors concerning the fragmen- 
tation progress, this has to be done manually. We enhanced the corresponding 
routines to report some informations such as the ^f, which enable us to cal- 
culate each string breaking point and all meeting points in four dimensions. 
Please recall, that this is all the information we need. We have checked, that 
the proper times reconstructed from these values match those used in the 
fragmentation process. 

As an example we consider the fragmentation of a g^fg-string with Jetset. 
For simplification purposes we switched off transverse momentum production 
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Fig. 2. Distribution of production (upper panel) and formation points (lower panel) 
for the example Pq^q = (±10, 0) GeV, Pg = (0, 5) GeV (left) and Pg = (3, 4) GeV 
(right). Dark regions refer to larger decay and formation probabilities. 

at the string breaking vertices. On the left hand side of Fig. 2 we show the 
spatial distributions of the string breakings, i.e. production points, and the 
yoyo-formation points for a system where the quark and antiquark momenta 
(±10 GeV) are along the x-axis (with opposite directions) and the gluon moves 
perpendicular to them into the ^/-direction (5 GeV). The spatial boundaries of 
the production points are set by the initial (anti-) quark and gluon momenta. 
Since the gluon looses energy to both the adjacent string pieces the production 
points have transverse coordinates y < 2.5 fm. As one might have expected, 
the yoyo-formation points are slightly shifted outwards compared to the pro- 
duction points. Taking the same parameters, but slightly tilting the gluon 
momentum out of the transversal direction, the spatial distributions change 
according to the right hand side of Fig. 2. Please note however, that these are 
more or less pathological examples, since in real event simulations (some hard 
scattering with initial and final state radiation) string regions with invariant 
masses as large as chosen here are quite rare. 

The sharp structures, best seen in the distribution of the production points for 
the tilted gluon (Fig. 2, right top), are due to mass constraints and only visible, 
if one neglects transversal momenta, as done in this example: Structures caused 
by different hadron masses m are smeared if one has to switch to the transverse 
mass mj_. We point out, that we have excluded production and formation 
points, which are due to a faulty reconstruction caused by the abovementioned 
simplifications in the Jetset code. These would give rise to spurious points at 
transversal components larger than what is allowed by the gluon component 
and contribute with less than two percent to the total yield. 



In both cases one clearly identifies a new string region as the horizontal dis- 
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Fig. 3. Production (top) and formation (bottom) points for a typical Hermes event: 
(u = 14GeV, = 2.5 GeV^). The target nucleon is located at the origin, the virtual 
photon is coming in from the left. Labels in the plots correspond to explanations in 
the text. 

tribution displaced from the This region appears in the time-evolution 

of the string after the gluon has lost all of its momentum to the string. 

We now apply our method of extracting four-dimensional vertices to two prac- 
tical examples. The first is an electron (virtual gamma) induced reaction on 
a nucleon at HERMES energies and the second one is a pp-coUision at RHIC 
energies. 

The experimental data on 7*p and 7*^4 reactions at HERMES imply some 
complicated kinematic cuts, e.g. in the zz-Q^-plane where u denotes the pho- 
ton's energy and its virtuality. In this study we will neglect this com- 
plication and simply replace the whole averaging over a multi dimensional 
parameter space by using some mean values. Since the underlying physical 
processes vary drastically (e.g. importance of resolved photon-nucleon inter- 
actions) within the full experimentally accessible kinematical region, our ex- 
ample is not representative for a// possible HERMES events. In Fig. 3 we show 
the spatial distribution of production and formation points for a 7*p-reaction 
at (z/) = 14GeV and (Q^) = 2.5 GeV^ This corresponds to an invariant mass 
(W) — 5 GeV of the photon-nucleon system. The photon is coming in from the 
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left and strikes a nucleon which is located at the origin. While PYTHIA de- 
scribes all spectra of non-charmed hadrons at HERMES energies astonishingly 
well (jsl: IgI: lllh . its ability to describe charm production has never been tested 
in this kinematic regime. We therefore exclude all charmed hadrons from our 
analysis. Note that both production and formation points lie mainly inside a 
radius of ~ 6fm around the interaction point. If one compares these values 
with a typical diameter of a complex nucleus one might expect strong (pre- 
)hadronic FSI after the (production) formation time in case of photonuclear 
reactions. 

In contrast to common belief the DIS does not simply lead to a one- dimensional 
string configuration. Instead of that, the distributions in Fig. 3 show some 
(elliptical) structures which are again due to non vanishing mass terms and 
need some more detailed explanations: Due to kinematics all production and 
formations points lie within circles of radius W/2 in the cm-frame of the 
photon-nucleon system. By boosting to the lab frame one obtains the observed 
ellipses with maximum transverse dimension W/2 and extremal longitudinal 
values of x~ ~ — 0.4fm and ~ 14.5 fm. 

The transversal distribution of the vertices is caused by two mechanisms. The 
first one is the transverse momentum p± generated in the string breakings 
that leads to a small transverse spread which is best seen at longitudinal 
distances around 7fm. The by far more important effect, however, stems from 
the intrinsic transverse momentum k± of the nucleon constituents. This can 
easily be understood by considering the simple scenario, where the virtual 
photon gets absorbed by a single quark inside the proton. In that case a 
single string is spanned between a diquark at the origin and the struck quark. 
Without intrinsic momenta, such a string would simply expand along the 
initial photon direction. The intrinsic kj_ rotates the direction of the string 
(without changing its invariant mass) and leads to the transversal distribution 
of production and formation points. (We used the default value (fc^) = 1 GeV^ 
of Pythia in our simulations.) 

The yoyo-formation point of the first-rank hadron, i.e. the hadron labeled 
'1' in Fig. 1, is given by + ■ the case of one single quark-diquark 
string its longitudinal coordinate in the lab frame is therefore confined to 
the narrow region between 14.1 fm and 14.5 fm if one neglects all transverse 
momenta. Correspondingly, the formation point of the highest-rank hadron, 
i.e. the hadron labeled '3' in Fig. 1, is given by x~ + ^n^~^ and, hence, its 
longitudinal coordinate lies between — 0.4fm and 14.1 fm. The finite invariant 
mass W of the string gives rise to a lower boundary [C,^ > [mh/WY) for 
the possible ^-values. This lower limit leads to the two distinct branches that 
are visible in Fig. 3: While the lower limit is negligible for pions, it becomes 
important for the heavier mesons and baryons; as a consequence the latter 
have formation points with larger (here: positive) longitudinal component. A 
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Fig. 4. Distributions of the production times ti, t2 and the formation time 
tp for a virtual event. The kinematics are the same as for Fig. 3 

{u = 14GeV, Q2 = 2.5 GeV^). 

similar (but more involved) mechanism is at work for the production points 
and leads to the strong (dark) band at low longitudinal positions. 

The importance of the production points for hadron attenuation in nuclear 
reactions has been extensively discussed in Refs. (0; 0; S H 0)- The con- 
stituents of a hadron that is formed at point Fi in Fig. 1 are produced at 
the string-breaking vertices Vi and Vi_i. We denote the larger of the two pro- 
duction times ^2 and the smaller one ti. At the later production point, i.e. at 
time t2, the string splits into a color neutral object (prehadron) - which has 
the four-momentum of the final hadron - and a remainder string. The final 
state interactions of this prehadron with the nuclear environment modify the 
spectrum of the ultimately detected hadrons but do not influence the actual 
string fragmentation process. Fig. 4 shows the distributions of the produc- 
tion times ti and ^2 of the first and second hadron constituents as well as 
the yoyo- formation times tp- The structures in the time distributions simply 
reflect the complex features that already showed up in Fig. 3. The presented 
distributions might be used as input for various attenuation studies such as 
the recent prehadron absorption model of Ref. |l7[). 

In Fig. 5 we show the average production times (ti) and (t2) as well as the 
average formation time (tp) as function of the fractional hadron energy Zh = 
Efjv. One clearly sees that the hadron formation times increase with energy 
whereas the production times of prehadrons become small if they carry a large 
energy fraction Zh. According to the Lund model the prehadronic final state 
interactions therefore dominate the attenuation of high energy hadrons in DIS 
of nuclei. 

For comparison Fig. 5 also shows the simple time estimates of Ref. (jlil ) that 
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Fig. 5. The solid curves show the average production and formation times of hadrons 
in certain bins of the fractional hadron energy Zh = Eh/v. The kinematics are 
the same as for Fig. 3 {v = 14GeV, = 2.5 GeV^). The dashed curves repre- 
sent the results of the analytic estimate derived in Ref. for the decay of an 
one- dimensional string into infinitely many fragments. 

were derived analytically for the fragmentation of an one- dimensional string 
into infinitely many fragments. Obviously, the analytic estimates yield only a 
rough approximation of the average production and formation times in Jet- 
set. Furthermore, it is not enough to know the average production and forma- 
tion times if one wants to perform realistic calculations of nuclear reactions. 
Our method allows for a complete event-by-event simulation assigning the 
production and formation points and times to all hadrons in each Jetset 
event. 

The distribution of production and formation points for pp-collisions at RHIC 
(a/s = 200 GeV) shown in Fig. 6 can be easily understood: Around the interac- 
tion point - corresponding to mid-rapidity - one has the broadest transverse 
momentum distribution and hence the largest transverse coordinates of the 
production and formation points. The regions with rapidity y ^ oo mainly 
consist of remnants, that lead to hadron production along the beam axis. 
Since the rate of particles dN/dp±dy decreases approximately exponentially 
with increasing pj_, also the distribution of production and formation points 
in transversal direction does. (This is in contrast to the HERMES regime, 
where the photon sets an upper limit on the energy of the particles. Here 
this upper limit far out of reach and the transverse distribution is practically 
unconstrained.) 

Of major experimental interest are high p^-hadrons (e.g. p± > 4: GeV) at 
mid-rapidity {\y\ < 0.5). The production and formation points applying these 
cuts are shown in Fig. 7. While the distribution of formation points supports 
the popular belief that high-p^ mid-rapidity hadrons at RHIC are produced 
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Fig. 6. Production points for RHIC: N + N@^/s = 200 GeV. The picture is sym- 
metric in horizontal direction, the contribution extends up to ±100 fm. 

outside any fireball or deconfined phase, the distribution of production points 
tells another story: Most of the production points are within a transversal 
distance below 2fm. However, we point out that the high hadron densities 
created in ultra-relativistic heavy ion collisions might strongly influence the 
string fragmentation process and that in the presence of a deconfined phase 
(QGP) the whole concept of hadronic strings becomes meaningless. For heavy 
ion collisions the presented production times and lengths should therefore be 
rather seen as lower estimates. 



In order to quantify the findings from above. Tab. 1 presents the average 
proper times (r) of the production and formation vertices as well as the average 
production and formation times (t) in the laboratory frame. For comparison 
we also list the average values of the fraction t/r and the Lorentz boost factor 
7 = Eh/mh of the considered hadron. Note that for each hadron there exist two 
production (proper-)times (ri) ti and {T2) t2 - corresponding to the production 
of the first and second hadron constituent - as well as one yoyo-formation 
(proper-)time [tf) tp- The larger of the two production times t2 is also called 
prehadron production time (UI^. 

The average production and formation proper times at HERMES seem to be 
almost independent of the particle species. On the other hand, one observes 
the same universal feature for the average laboratory times at RHIC conditions 
(with the implemented cuts). Furthermore, the ratio (t/r) deviates from the 
Lorentz factor (7) of the corresponding hadron. Hence, the times that we have 
extracted from Jetset are in contradiction to the widely-used assumption 
that the hadron production and formation times in the lab frame are simply 
given as a product of the Lorentz-boost factor 7 and an universal constant 
proper time. 
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Fig. 7. Production and formation points for RHIC(N + N@y/s = 200GeVj for 
mid-rapidity (\y\ < 0.5) and high pj_: p± > 4GeV. 

In the presented work we have extracted the production and formation vertices 
of hadrons from the Lund fragmentation as modeled by Jetset. In contrast to 
other approaches, such as analytic derivations, we are able to extract the full 
four- dimensional information for each hadron produced in the fragmentation 
of a (gluonic) string on an event-by-event basis. For the first time we have pre- 
sented quantitative results for realistic photon-nucleon reactions at HERMES 
energies and pp-coUisions at RHIC using the Monte Carlo generator Pythia. 
However, we stress that our method can also be used for all other kinds of 
event generators which are based on Jetset for the fragmentation. We find 
that a large fraction of production and formation points at HERMES ener- 
gies lie inside typical nuclear radii which underlines the importance of (pre-) 
hadronic FSI. Also the production points of high-p^ mid-rapidity hadrons at 
RHIC seem to lie inside the fireball region. Furthermore, we conclude that the 
conventional determination of hadron production and formation times by sim- 
ply multiplying a (universal) constant proper time with the corresponding 7 
factor leads to misleading estimates for the hadron production and formation 
points. 



This work is thought to be a Jproof of concept'. Detailed transport-theoretical 
studies of nuclear reactions that explicitly model the (pre-)hadronic FSI 
in the nuclear environment will follow. We point out that any event-by-event 
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5.7 


13.8 
6.9 
6.0 


18.7 

8.5 

6.4 


51.04 
9.31 
7.12 



Table 1 

Average production and formation times for pions, p mesons and baryons within 
specified kinematics ("HERMES", "RHIC"). Proper times are labeled by t, while t 
represents the corresponding time in the laboratory frame. 7 is the relativistic boost 
factor of the hadron. 



simulation based on our findings can account for all sorts of experimental cuts 
and detector efficiencies. Such models can therefore be directly compared with 
experimental data or used for detector simulations. They will help to clarify to 
what extend the experimentally observed hadron attenuation in nuclear DIS 
and ultra-relativistic heavy ion collisions is caused by the rescattering of the 
produced colorneutral (pre-)hadrons after fragmentation and what is caused 
by additional in-medium phenomena like a modification of the fragmentation 
function, a partonic energy loss or recombination effects. 
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